Interchain Coupling Effects and Solitons in CuGeOa 
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The effects of interchain couphng on solitons and sohton lattice structures in CuGeOs are explored. 
It is shown that interchain coupling substantially increases the soliton width and changes the soliton 
lattice structures in the incommensurate phase. It is proposed that the experimentally observed large 
soliton width in CuGeOa is mainly due to interchain coupling effects. 
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The inorganic spin-PcijCrls (SP) system, CuGeOs, has 
attracted much attention^ recently. Pure CuGeOs has a 
SP transition at Tgp ~ 14.3K. Below Tgp, the system is in 
a dimerized spin singlet state, and the gap to spiii, triplet 
excitations is A ~ 24.5-KBB. It has been arguedflfl that 
the spin susceptibility measurements in this system can 
be understood only if the nearest neighbor Heisenberg 
Hamiltonian is augmented by a substantial next-nearest- 
neighbor (NNN) couphng, J2 ^ 0.2J, where J is the 
nearest neighbor coupling. The NNN coupling has been 
shown to be important in other experimentsQ as well. Re- 
cently, the soliton lattice structure in CuGeQ^-was stud- 
ied using neutron and NMR measurementaO'u. An im- 
portant observation is that the soliton width ^ ^ 13.6a, 
where a is the lattice-.spacing, is much larger than the 
theoretical predictionB for the unfrustrated SP chain, 
^0 = Ja7r/(2A) - 8a. 

The soliton width in CuGeOs can be changed by in- 
terchain coupling and frustration J2. Without inter- 
chain coupling, the soliton width is generally ^ oc v/A, 
where v is the spin wave velocity and A is the trijdei, 
gap in the SP chain. From nonlinear simga modeM, 
V = ^;n^/l — AJ-)/Ji . and the dimensionless coupling 



constantliil5s(J2) = 5s(0)/-\/l — 4J2/>/i- So, the frustra- 
tion will decrease the soliton width. One can show that 
the soliton width ^ = Jav/A. The triplet gap A will de- 
pend on both SP coupling and frustration. However, to 
compare with experiments, the triplet gap A ^ 24:.5K, 
while V decreases with increase of frustration, so the frus- 
tration will make the soliton width in a CuGeOa chain 
smaller than ~ 8a. Thus, the observed soliton width 
C/Co ^ 1.7 is not an effect due to the frustration J2 in 
CuGeOa. 

It is knjQwn that the interchain coupling in CuGeOa is 
not smallH: Jc ~ 120K, Jb ~ O.lJc, and Ja ^ -O.OlJc- 
The phase diagram of the unfrustrated SP system in the 
presence of inteijchain coupling was studied by Inagaki 
and Fukuyamalla. In the present paper we calculate the 
effect of interchain coupling on soliton width. Although 
the singlet SP ground state is destroyed by large inter- 
chain coupling, in the strongly dimerized regime weak 
interchain coupling cannot affect the spin singlet ground 



state. Similarly, weak interchain coupling has very little 
effect on the susceptibility above Tgp- From mean field 
theory, it can be shown that the effects are typically of 
the order of 3%. 

The low energy excitations involving broken spin sin- 
glets are more sensitive to interchain coupling. Once a 
singlet pair is broken in one chain, it is energetically less 
costly to break another pair in a neighboring chain due 
to the interchain coupling. Thus, the interchain coupling 
will affect the low energy excitations in the dimerized 
phase. For the same reason, the interchain coupling can 
lead to an increased soliton width: the local magnetiza- 
tion of a soliton in SP chain is M{x) ~ cosh~"^(a;/^) (c/. 
below), and the gain in the interchain coupling energy 
increases with the increased soliton width ^, hence an 
increased soliton width. 

In this paper, the effects due to the interchain cou- 
pling are studied using bosonization. We consider the 
Hamiltonian: 
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where Si,j (uij) is the spin operator (lattice distortion) 
of site-i in chain-j, and K is the spring constant. (3 is the 
spin-lattice coupling constant, ao = J2/ J '^s the strength 
of frustration and 7 = J±/ J is the interchain coupling 
constant. Note that we have used the adiabatic approxi- 
mation for the phonons. 

The interchaixi_coupling is treated in the mean-field 
approximationll3llj: S^j ■ Si^j+f, = Sfj{Slj^^). Thus, 
the Hamiltonian is transformed into a single chain 
model and the interchain coupling term transforms into 
-Y.,lJZh,St + llJZY^^hl where h^ = {Sl^_,^) = 

— {Z/Z){S^j) is the effective field. We have introduced 
the notion of an effective coordination number Z to 
account for {Sfj_^_^) ^ —{Sfj). In the soliton lattice 

phase, Z is site dependent. The effective Hamiltonian is 
then transformed to a fermionic Hamiltonian using the 



Jordan- Wigner transformation, and the fermionic model 
is bosonized. The resulting continuum Hamiltonian islU: 
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--^f dxh{x) cos{2(f>) + I dx^h^ix), (2) 

where a is the short distance cutoff, u{x) = (— l)*Wi/a 
and h{x) = {—lyhi/a are slowly varying variables. The 
bare value of the Umklapp term is 53 = 1 — 3ao and we 
shall assume that the renormalized 53 has approximately 
the same ao'dependence. 

Here v is the spin wave velocity and Kp is a crit- 
ical exponent. For the unfrustrated AFM Heisenberg 
spinj-chain, the following renormalized values should be 
usedllj: v = tt/2 and Kp = 1/2. However, the values 
of V and Kp in the frustrated Heisenberg model is de- 
pendent on ao. With increase of frustration, the spin 
wave velocity decreases. In nonlinear sigma modeO, 

V = wqVI ~ 4q;o for small ao. In a recent numerical cal- 
culation, it is shown that the variation of renormalized 

V with ao is approximatelyllj p^ = 2v/7t ~ (1 — 1.12ao). 
Although we don't know Kp at finite ao precisely, it is 
knowiO that Kp > 1/2 at small frustration ao < ag ~ 
0.3. Define pK = 2Kp — 1, then px -^ 1 and pv < 1. 

For the ground state, u{x) = uq, h(x) — ho, and Z = 
Z. If there is no Umklapp interaction 33, the model (g) is 
exactly soluble by Bethe AnsatzES. Due to the (73-term, 
an approximation has to be employed. We use the self- 
consistent Gaussian approximation, which was previously 
shown to be reliablecl for the SP chains. To this end, we 
write (t){x) = 0s (a;) + 0(a;), where ^^ is the semiclassical 
solution and (j) is the fluctuation around </)<,. We retain 
(j) to quadratic order in the action. The first order term 
vanishes at the saddle point. Then, 



92(40,) + — sin(40. 



where 



and 



e 






^__0^ 2_853 4 
£,1 a^K "" a^" - 







The quantity a is given by 



(T = e ^^ ' . 



(3) 
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a describes the renormalization of 51 and 53 due to fluctu- 
ations of (j) around </), . Note the derivation is appropriate 



for systems with a finite gap, otherwise, a has infrared 
divergence. In deriving these equations, we have used the 
self-consistent equations 



u{x) = {(3J/K)asin{2(j),), 
h{x) = (Z/Z)crcos(20s)- 



(7) 
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The "uniform" ground state (i.e. Ui — (— l)'uoa and 
hi — (— l)*/ioa) calculated from Eq. (g) corresponds to 
sin(4(/)s) = 0. Therefore, either 0s = or 7r/4. When 
0s = 0, u{x) = 0, ho ^ 0, the system has long range 
AFM order; when 0s = 7r/4, h{x) = 0, uq ^ 0, the 
system is dimerized. The spin-lattice coupling and inter- 
chain coupling strengths dictate the solution that has the 
lowest energy. It is clear that on symmetry grounds these 
two homogeneous solutions are mutually exclusive. 

The calculation of ground state energy and excitation 
spectrum in the self-consistent Gaussian approximation 
is straightforward and can be read off from RefJj, so we 
will just state the results. The excitation spectrum cor- 
responding to uniform 0s is 



u!q ~ Jav\l q-^ 
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In the dimerized phase 0s — 7r/4 and 0o — 0, and it can 
be easily shown that go — l/^o- From Eq. (||), the spin 
triplet excitation gap is A = Javqo. The scaling of the 
spin triplet excitation gap is 
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where S = (3u is the bond alternation induced by lattice 



dimerization. For the unfrustrated SP chain, K^ 



1/2, 



Pk = 0, and we obtain the correct scaling relatior£3 A c>c 
,52/3. 

Comparing the ground state energies between dimer- 
ized phase (0s = 7r/4) and AFM phase (0s = 0), we get, 
for Pk ^ 1, that the crossover from 0s = to 0s — 7r/4 
is determined by 



Ck 



_ JjZ / A \PK 
Apv VTr^Jp, 



= 1. 
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If Ck < 1 then 0s = 7r/4, 6 = /3u, h — 0, and the sys- 
tem is dimerized; otherwise, S = "fZh, u — 0, and the 
ground state has long range AFM order. A similar crit- 
ical value of the interchain coupling for the unfrustrated 
SP model was derived by Inagaki and Fukuyamali3. Us- 
ing A ~ 24.5i^ and J ~ 120i^, we find that Ck ^ 
0.98pv{0-021/ pv)''^ < 1, which is consistent with the fact 
that CuGe03 is dimerized at zero temperature. Since 
Pv (pk), which is less (greater) than 1(0), decreases 
(increases) with increase of frustration, the frustration 
pushes the SP-Neel transition to stronger interchain cou- 
pling. If pi, ^ 1 and Pk <^ 1, the interchain coupling is 
close to the critical value for the long range AFM order 
to be realized. This is probably significant, as materials 
with only 3% Zn exhibit long range AFM order. 



Consider the nonuniform solutions of Eq. (g), which 
correspond to the soUtonic excitations of the dimerized 
phase. In zero magnetic field, the soliton solution has 
lower energ5cl than the spin-triplet excitation with the 
gap A. In this case, we can continue to assume that Z jZ 
is a constant, and the soliton solution from Eq. (J^) is 



interchain coupling is to renormalize ^, as in Eq. ( [T^ ) 
The periodic solution of Eq. 



u(x) (X sin(2(/)s) = ±tanh 
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with the half-width 
Jav 
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The corresponding magnetization, M[x), is 

M(.t) =cosh"i(a;/C)/(27rO (14) 

and the staggered magnetization, N{x), is 



Nix)^{-rJ^cosh-\x/0. 
V Ztt 



(15) 



From Eq. (|l^), we can see that the soliton width is in- 
creased by interchain coupling, while the effects of frus- 
tration ao — J2/J appear only in A and v. 

Without interchain coupling, the soliton width is ^0 = 

JaTTpJ{2A) < 8a for CuGcOa. With Z = Z, {Sl^^ = 

'{Sfj), the soliton excitation described by Eq. (12) is 



an array of antifcrromagnctically coupled solitons in a di- 
rection perpendicular to the SP chains. As discussed in 
the introduction, the spins are antiferromagnetically or- 
dered inside the soliton, therefore increased soliton size 
will lead to a gain in the interchain coupling energy. This 
is the physical origin of the increased soliton width in 
the presence of interchain coupling. When the interchain 
couphng is so strong that the soliton width diverges, i.e. 
Ck — 1, the system crosses over to the long range anti- 
ferromagnetically ordered phase. This criterion Ck = 1, 
derived from consideration of the soliton size, is the same 
as that derived from Eq. ( [ll| ) for small px- 

If the solitons are excited locally without forming a 
regular array, {S^j_^_^) 7^ —{S^j), and Z/Z is smaller 
than unity. The soliton size will still be increased with 
the increased interchain coupling 7, and AFM domains 
will form and increase in size. We propose that for the 
solitonic excitations due to magnetic or non-magnetic im- 
purities similar AFM domains will form. The percolation 
of these AFM domains must be relevant to impurite-|iB|- 
duced AFM ordering in the doped CuGeOs systemsc2lEil. 

In a high magnetic field, there will be a transition 
from the zero field commensurate dimerized phase to 
the incommensurate soliton lattice phaseEd. The ground 
state of this soliton lattice phase corresponds to the self- 
consistent periodic solutions of the coupled nonlinear 
equations similar to Eq. (0). Here, we consider only 
the qualitative aspects of the effect of the interchain cou- 
pling in the soliton lattice phase. To this end, we assume 
a constant ^ in Eq. (H) for each chain. The effect of the 



is: 
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where ^ is the soliton width defined in Eq.|-(tf3[), and sn, 
dn, and en are the Jacobi elliptic functionsErTThe inter- 
soliton distance is 2K{k)k£^ [see Fig.(|l|)], where the mod- 
ulus k of the elliptic integral K is related to the total mag- 
netization induced by the external magnetic field-that can 
be derived from the minimization of the energyE3. These 
equations are the same as,those of the soliton lattice with- 
out interchain couplingEil. Close to the commensurate- 
incommensurate transition, M{x) ^ ^ |iV(x)|, we 
can assume that these equations still approximately de- 
scribe the soliton lattice in the presence of interchain 
coupling: each chain has a magnetization M(x)^ and 
the neighboring chains have staggered magnetizations 
N{x) and —N{x). The effective coordination number 
Z/Z is smaller than unity. Thus, Z/Z = constant is 
a reasonable approximation close to the commensurate- 
incommensurate phase transition, where most of the ex- 
periments are carried out. |-. 

In the experiment of Kiryukhin et alx\, the value of ^ 
is measured in the soliton lattice phase. The measured 
value of ^ is ^ = (13.6 ± 0.3)a, which is substantially 
larger than the value ^0 ~ 8a without interchain cou- 
pling. If we use an average Z/Z, from ^ — (13.6 ± 0.3)a 
we estimate Z/Z ^ 0.8 at ao = 0. With finite aoi 
the estimated Z/Z will be larger. In these experiments 
the inter-soliton distance is 40 — 70a, so fc ~ 0.7 — 0.9. 
The magnetization M{x) and staggered magnetization 
N{x) are shown in Fig. (|l]) for k = 0.8. The ratio 
{M ix)) / \J {N"^ {x)) ~ 0.1, which is the right order of mag- 
nitude to giwe effectived Z/Z - 0.8 - 1.0. In the NMR 
experimental, the maximal value of the effective local 
spin is measured to be Smax ~ 0.065, which is smaller 
than the value corresponding to a single chain given by 
^A/(j7r2p^) > 0.14. However, because the solitons in 
the neighboring chains are antiferromagnetically coupled, 
the effective local magnetization at each site is reduced 
by the interchain effects; a reduction factor of the order 
of 2 — 3 of S'inax is quitc reasonable. 

In conclusion, we have studied the combined effects of 
the frustration and interchain coupling on solitons and 
soliton-lattice excitations in CuGeOa. We find that the 
interchain coupling can substantially increase the soliton 
size, while the NNN frustration will decrease it. When 
the interchain coupling strength is close to the SP-Neel 
transition, the size of the soliton diverges. The analysis 
of the soliton structure prcsciited here is consistent with 
the experimental observationaJ^B. 
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FIG. 1. The lattice distortion and magnetizations in the 
soliton lattice phase with f = 13.6, k = 0.8. The scale of u{x) 
is arbitrary. 
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